Gravitational wave generation in loop quantum cosmology 
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We calculate the full spectrum, as observed today, of the cosmological gravitational waves gen- 
erated within a model based on loop quantum cosmology. It is assumed that the universe, after 
the transition to the classical regime, undergoes a period of inflation driven by a scalar field with a 
chaotic-type potential. Our analysis shows that, for certain conditions, loop quantum effects leave a 
clear signature on the spectrum, namely, an over-production of low-frequency gravitational waves. 
One of the aims of our work is to show that loop quantum cosmology models can be tested and 
that, more generally, pre-inflationary physical processes, contrary to what is usually assumed, leave 
their imprint in those spectra and can also be tested. 
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I. INTRODUCTION 

Although gravitational waves of cosmological origin 
have not yet been detected, they are at present the object 
of a considerable research effort, as they may provide us 
with a unique telescope to the earliest stages of the for- 
mation of the universe. At the same time we also witness 
an increasing interest in the applications of the ideas of 
loop quantum gravity (for a review, see Ref. [l|) to the 
problems of cosmology, a field known as loop quantum 
cosmology, interest that started after a series of seminal 
papers by Bojowald where a number of important 
results were obtained, among them the possibility of re- 
moving in a natural way the presence of the cosmological 
singularity (for a review about loop quantum cosmology, 
see Ref. [|). 

Tests of l oop quantum cosmology have already been 
proposed [4|-[lO|. showing that loop effects may appear, 
albeit in an indirect way, on the cosmic microwave back- 
ground radiation on the largest scales. Loop quantum 
cosmology gives rise to changes in the dynamical equa- 
tions driving the expansion of the universe; these are con- 
nected with modifications in the equation of state of the 
matter content of the universe, which in turn result in the 
production of gravitational waves. These gravitational 
waves are the focus of our work, where we show that they 
may leave an important imprint in today's power spec- 
trum. This at first may seem surprising, as it has usu- 
ally been assumed that inflation, by its characteristics, 
among them the enormous increase in the scale of the 
universe, would remove any kind of information coming 
from physical phenomena taking place in pre-inflationary 



times. This is not the case. Pre-inflationary physical fea- 
tures affect in a different way different frequencies and 
the memory of such differences survives through the in- 
flationary period, and can be read today in the power 
spectrum. 

Our paper is organized as follows. In the next sec- 
tion we describe the loop quantum cosmology model used 
in our work and write the equations of motion for the 
semiclassical and classical stages of evolution of the uni- 
verse. Taking into account the constraints arising from 
measurements of the cosmic microwave background ra- 
diation, we specify the values of the parameters and the 
initial conditions and solve numerically the evolutionary 
equations, from the semiclassical prc-inflationary epoch 
till the present time. In Section IIIII we calculate the 
full gravitational- wave spectrum, for frequencies ranging 
from about 10^^'' rad/s to about 10^ rad/s, using the 
method of continuous Bogoliubov coefficients. The influ- 
ence on the spectrum of the ambiguity parameters of loop 
quantum cosmology is carefully analyzed. In section IIVI 
we summarize the results we have obtained. 

II. THE EVOLUTION OF THE UNIVERSE 

The evolution of the universe is naturally divided in 
two stages, corresponding to the semiclassical and the 
classical regimes. During the semiclassical stage of evo- 
lution, the equations of standard cosmology have to be 
modified in order to account for loop quantum effects. 
After a few Planck times of evolution, the transition be- 
tween the semiclassical and classical regimes takes place 
and the universe further evolves according to the stan- 
dard cosmological model, namely, undergoes a period of 
inflation, driven by a scalar field (j) with chaotic-type po- 
tential 
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followed by reheating and, successively, by radiation- 
dominated, matter-dominated and dark energy- 
dominated periods. 

During the semiclassical stage, the modified evolution- 
ary equations for the scale factor a and the scalar field cj) 
are given by 0, [III, 
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where a flat Friedmann-Robertson- Walker background is 
assumed, mp is the Planck mass, a dot denotes a deriva- 
tive with respect to the cosmic time t and the following 
notation was introduced llSlI 
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FIG. 1; Time evolution of the scalar field <j> for the case j = 
100 and i = 3/4. Because of loop quantum effects, the scalar 
field increases naturally from its initial value of about 3.6 x 
10^^ mp to a maximum value of about 3 mp , which guarantees 
a long enough inflationary period. 
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In the above expressions I and j are the so-called ambi- 
guity parameters and denotes the Planck length. 

Equations ^ and ^ are solved numerically for the 
following values of the parameters and the initial con- 
ditions: = 10~^mp, Oi = ■y/7£p, = 3.6 x 
10^^j~^^/^mp and 0i = 2 x 10^*' mp, with being fixed 
by Eq. ([2]) . This equation is also used to check the accu- 
racy of the numerical solution. The initial value for scalar 
field was chosen such that the uncertainty principle Q , 
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^ Different values for the Barbero-Immirzi parameter 7 can be 
found in literature. We use the value obtained by Meissner from 
black-hole entropy calculations Il4ll . 



is marginally satisfied. 

The ambiguity parameter i is quantized, taking values 
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The other ambiguity parameter, j, takes half- integer val- 
ues greater than one. However, if one demands the initial 
value of the scalar field, to be much smaller than the 
Planck mass (say, < 10"'^ mp), then the parameter 
J is constrained to be much bigger than one, namely, 
j > 87.3. In what follows, we will consider j to be a 
continuously varying parameter with values greater or of 
the order of 100. 

As our numerical calculations show, after a short pe- 
riod of time, the functions Di^q) e F{(i) approach their 
classical values, namely, _D = 1 and F = 0, and the 
semiclassical corrections can be neglected in Eqs. 
(|4]). The resulting simplified evolutionary equations are 
then solved up to the end of the infiationary period. 

At the end of the semiclassical stage of evolution, the 
scalar field increases from a small value (much smaller 
than the Planck mass) to a value of the order of one 
Planck mass (see Fig. [1]). This is due to the fact that 
the first term on the right-hand-side of Eq. (|4]) acts as an 
anti-friction term, pushing the scalar field up the poten- 
tial. During the subsequent classical stage of evolution, 
the scalar field continues to increase for a while, reach- 
ing a maximum value of about 3 mp (for j = 100 and 
£ = 3/4). This value of ^max is enough for the universe 
to expand about 60 e-folds during the inflationary period. 
Therefore, loop quantum effects can set the initial con- 
ditions for successful chaotic inflation in a natural way 

At the end of the inflationary period, the scalar field 
begins to oscillate around the minimum of the potential, 
transferring its energy to a radiation fluid, thus reheating 
the universe. The decay of the scalar held into radiation 
is governed by a dissipative coupling introduced in the 
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evolutionary equations, which now read 
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where pr is the energy density of radiation and is the 
dissipative coefficient. Since any preexisting radiation 
fluid would have been diluted during the inflationary pe- 
riod, we choose the energy density of radiation at the 
beginning of reheating to be zero, = 0. For the dis- 
sipative coefficient we choose = 2 x 10~^mp. 

After a while, the energy density of the scalar field be- 
comes much smaller than the energy density of radiation, 
meaning that the former can be consistently neglected. 
The evolutionary equations then become^ 
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where ao is today's value of the scale factor and pr.o 
4.13x10-" J/in, p^^o = 2.34x10-1" J /in and pde,o = 
6.20 X 10~^'^ J/m'^ are, respectively, today's values of the 
energy density of radiation, matter, and dark energy. For 
these values of p,-^o, Pm.o, and pde,o> the value of the Hub- 
ble constant is Hq = 71 km s~"'^Mpc~"'^. We take the 
value w = —I for the equation-of-state parameter of dark 
energy. 

To finish this section, let us point out that for a given 
value of the requirement that the universe expands 
enough during the infiationary period (at least 60 e- 
folds), imposes a lower bound on the value of £ (see 
Fig. [2]). 



III. THE GRAVITATIONAL- WAVE SPECTRUM 

Gravitational waves are generated in an expanding uni- 
verse, giving rise to a spectrum extending over a wide 



^ In the previous stages of evolution we have used the natural 
system of units, with h = c = 1 and mp = G~^^^ = 1.22 X 10^^ 
GeV, while here we are using the international system of units. 
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FIG. 2: The shaded region corresponds to values of the am- 
biguity parameters £ and j for which the universe expands at 
least 60 e-folds during the infiationary period. Such an ex- 
pansion is achieved if the scalar field grows from <j)i <^ mp 
to 0max ^ 3mp during the pre- inflationary epoch. The 
requirement that (fi ^ IQ-^mp, together with the choice 
= 10"^ mp and (pi = 2x 10"^ mp, imposes the constraint 
j > 87.3. 



range of frequencies |1 5142 If . In this section we calcu- 
late the full spectrum of the gravitational waves gener- 
ated within the loop quantum cosmology model described 
above. 

Loop quantum effects introduce modifications not only 
to the dynamical equations of evolution ©-([l]), but also 
to the equation for tensor modes. However, in this paper, 
the latter will be neglected, allowing for a considerable 
simplification of the calculations required to determine 
the energy spectrum of gravitational waves, while keeping 
unchanged the main features of the spectrum. 

The tensor perturbations hij to the Fricdmann- 
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can be expanded in terms of plane waves 



(27r)3/2a 



i(r;)V2fc 



K(r?,k)£,,(k,p)e*''-e(r7,fc)+H.c.] 



(17) 



where G is the gravitational constant, p runs over the 
two polarizations of the gravitational waves, fc = |k| = 
27ra/A = au is the co- moving wave number, ap is the 
annihilation operator, and is the polarization tensor. 
The mode function ^(77, k) obeys the equation of a para- 
metric oscillator. 
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where the potential U is given by 
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and a prime denotes a derivative with respect to confor- 
mal time ?/. 
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FIG. 3: Potential ?7 as a function of (cosmic) time t for the 
case I = 3/4 and j — 100. The inflationary barrier is located 
on the right (for t > W^tp), while the much smaller barrier 
arising due to loop quantum effects is located on the left (for 



For fc^ < [/, the above equation describes the produc- 
tion of gravitational waves with frequency cj = k/a, while 
for ^ U the equation is that of an harmonic oscillator, 
implying that no gravitational waves are produced. 

Within the standard model of cosmology, the poten- 
tial U has a pronounced barrier at the time of inflation, 
giving rise to a copious production of gravitational waves 
with frequencies up to the gigahertz. Within loop quan- 
tum cosmology, besides the above mentioned inflationary 
barrier, the potential U has another barrier at very early 
times (see Fig. [3]). The presence of this barrier leads to 
the creation of extra gravitational waves of low frequency. 
This barrier is slightly more pronounced, if one takes into 
account loop quantum corrections to the equation of ten- 
sor perturbations. 

In order to calculate the energy spectrum of the 
cosmological gravitational waves generated during the 
evolution of the universe, from the semiclassical pre- 
inflationary stage of evolution till the present time, we 
use the method of continuous Bogoliubov coefficients. 
This method, first applied by Parker [2^ to particle pro- 
duction in an expanding universe and then extended to 
the case of gravitons |23l - [25| . can be summarized as fol- 
lows (for applications of this method to several cosmo- 
logical models see Refs. (26l - [30j ). 

The gravitational- wave spectral energy density param- 
eter. Slew, is defined as 
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where h is the reduced Planck constant, G is the gravita- 
tional constant, c is the speed of light, H is the Hubble 
parameter, a; is the gravitational-wave frequency, /3 is a 
Bogoliubov coefficient and the subscript denotes quan- 
tities evaluated at the present time. 

The angular frequency loq takes values ranging from 
about 1.4 X 10~^^ rad/s (corresponding to a wavelength 
equal, today, to the Hubble distance) to about 10^ rad/s 



FIG. 4: Gravitational-wave spectra for £ = 3/4 and £ = 5/6 
(dashed line). In both cases j — 100. 



(corresponding to a wavelength equal to the Hubble dis- 
tance at the end of the inflationary period). 

The number of gravitons at a certain time is given by 
the squared Bogoliubov coefficient, 

|/3(i)p = i[x(t)-y(t)][x*(t)-r*(t)], (21) 

where * denotes complex conjugate and the functions X 
and Y are solutions of the system of differential equations 
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X, (23) 



which is integrated with initial conditions X{ti) = 
Y{ti) = 1. corresponding to the absence of gravitons at 
the beginning of the evolution. In the above system of dif- 
ferential equations, the scale factor, a{t), and its first and 
second derivatives, d{t) and a(t), are determined from the 
evolutionary equations presented in the previous section, 
namely, Eqs. ([2])-(|4]) for the semiclassical stage of evolu- 
tion and the infiationary period , ([TU | - p^ for the reheat- 
ing period, and (fT4)) -([T5 |) for the radiation-dominated, 
matter-dominated and dark energy-dominated periods. 

Using the above outlined method of continuous Bogoli- 
ubov coefficients, we can calculate the gravitational- wave 
spectrum for different values of the ambiguity parame- 
ters J and £. Let us first analyze the case of fixed j (say, 
j — 100) and varying £. 

For £ — 3/4, loop quantum effects leave a clear signa- 
ture on the spectrum, namely, an over-production of low- 
frequency gravitational waves (see Fig. |4]). For £ — 5/6 
the spectrum shows no influence of these effects. This can 
be understood as follows. Today's frequency of a gravi- 
tational wave that crossed the Hubble horizon at time t 
is given by 2nH{t)[a{t)/a{to)]- For i = 3/4, gravitational 
waves generated in the early universe {t < 10^ tp) have 
frequencies, today, of the order of (10~^^ — 10"^'') rad/s, 
corresponding to wavelengths of the order or smaller than 
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FIG. 5: Today's frequency of a gravitational wave as a func- 
tion of the time at whicii this wave first crossed the Hub- 
ble horizon. For £ — 3/4, gravitational waves generated 
in the early universe have frequencies, today, of the order 
of (10"" - 10"") rad/s (shaded region), corresponding to 
wavelengths of the order or smaller than the Hubble hori- 
zon today. For £ — 5/6, gravitational waves generated in 
the early universe have, today, frequencies of the order of 
(10~^® — 10~*^) rad/s, corresponding to wavelengths much 
bigger than the Hubble horizon today. For both cases today's 
maximum frequency of the gravitational waves is of the order 
of lO'' rad/s, corresponding to waves that crossed the Hubble 
horizon at the end of the inflationary period (t ~ 3 x 10^ ip). 



the Hubble horizon today (sec Fig. [S]). Therefore, these 
gravitational waves leave their imprint on the spectrum 
at low frequencies. For £ ~ 5/6, gravitational waves 
generated in the early universe have, today, frequencies 
of the order of (10~^^ — 10"'''^) rad/s, corresponding to 
wavelengths much bigger than the Hubble horizon today. 
Therefore, these gravitational waves leave no imprint on 
the spectrum"^. Note that in the case £ = 3/4 the scale 
factor grows about 60 e-folds during the inflationary pe- 
riod, while in the ease £ = 5/6 the growth of the scale 
factor is about 140 e-folds. 

For the case j = 2000 the situation is similar to the one 
described above (see Fig. [6]) . Loop quantum effects leave 
an imprint on the gravitational- wave spectrum only if the 
ambiguity parameter £ takes such a value {£ = 1/2) that 
the universe expands about 60 e-folds during the infla- 
tionary period. If the universe expands much more than 
60 e-folds (for instance, in the case £ = 3/4 the scale 
grows about 240 e-folds), the gravitational waves gen- 
erated before the inflationary period have wavelengths, 
today, much bigger than the Hubble horizon and, conse- 
quently, they leave no imprint on the gravitational-wave 
energy spectrum. 

We have been assuming that the ambiguity parame- 



^ The small rise on the spectrum at low frequencies for £ = 5/5 
is due to another effect, namely, an extra production of gravita- 
tional waves in recent epochs, when the universe became matter 
dominated and, then, dark-energy dominated. 
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FIG. 6: Gravitational-wave spectra for £ = 1/2 and £ = 3/4 
(dashed line). In both cases j — 2000. 



ter £ is quantized, taking values £ = 1 ~ l/(2n), with 
n S N. If, however, we consider ^ to be a free parameter, 
changing continuously from to 1, then our conclusions 
need to be slightly adapted. Namely, for each value of j 
there is a critical value of the ambiguity parameter, ^crit; 
for which the universe expands 60 e-folds during the in- 
flationary period, the minimum required in standard in- 
flationary cosmology. For values of £ close to ^crit, the 
gravitational waves generated during the pre- inflationary 
epoch leave an imprint on the spectrum at low frequen- 
cies. As we (continuously) increase £, this imprint shows 
up at lower and lower frequencies, completely disappear- 
ing when today's frequency of the generated waves is so 
low (ujQ < 10~^^ rad/s), that it corresponds to a wave- 
length greater than the Hubble horizon. On the con- 
trary, if we consider values of £ smaller than the critical 
one (in which case the scale factor does not grow enough 
during the inflationary period), the imprint of the pre- 
inflationary gravitational waves appears on the spectrum 
at higher frequencies. 

The above conclusions are illustrated in Fig. [T] For 
j = 100 the critical value of the ambiguity parameter £ 
is about 0.75 (see Fig. [3]). In this case, loop quantum 
effects leave their signature on the spectrum at frequen- 
cies ujo < 10"^^ rad/s. For £ = 0.70 and £ = 0.73 the 
gravitational waves generated prior to the inflationary 
period leave an imprint at frequencies higher than in the 
critical case, ujq < 10^^ rad/s and ojq < 10^^ rad/s, re- 
spectively. However, for such values of the ambiguity 
parameters, the scale factor grows about 45 and 54 e- 
folds, respectively, during the inflationary period, which 
is less than required by standard inflationary cosmology. 
For the £ = 0.77 (above ^crit = 0.75) the signature of loop 
quantum cosmology is located at frequencies lower than 
in the critical case, ujq < 10""'^^ rad/s. In this case, the 
scale factor grows about 71 e-folds during the inflationary 
period. 

To conclude, let us compare gravitational-wave spec- 
tra for a fixed value of £ (say, £ = 0.6) and varying j. 
For £ — 0.6; enough growth of the scale factor during the 
inflationary period is guaranteed for j > 560 (see Fig.[2|). 
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FIG. 7: Gravitational-wave spectra for I = 0.70, 0.73, 0.75 
and 0.77 {j = 100 in all cases). Loop quantum effects leave a 
signature on the spectra, namely, an over-production of gravi- 
tational waves at low frequencies. As I increases this signature 
is located at increasingly lower frequencies. For I = 0.77 the 
signature moved to frequencies so low, that it becomes almost 
unnoticeable in the spectrum. 
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FIG. 8: Gravitational-wave spectra for j — 100, 560 and 750 
{i = 0.6 in all cases). As j increases the loop quantum signa- 
ture shows up at increasingly lower frequencies. 



Therefore, as expected, for j — 560 loop quantum effects 
leave their imprint at low frequencies, loq < 10""'^^ rad/s 
(see Fig. [8]). For j ~ 750 the loop quantum signature 
moves to lower frequencies, ^ 10~^^ rad/s, while for 
j = 100 it moves to higher frequencies, coq < 10 rad/s. 
In the latter case, the imprint of loop quantum cosmol- 
ogy lies in the frequency band accessible to the Laser 
Interferometer Space Antenna (LISA) and, in principle, 
could be seen by this detector. Note, however, that the 
growth of the scale factor during the inflationary period 
in this case is just about 28 e-folds, which is manifestly 
insufficient within standard inflationary cosmology. 

To finish this section, let us point out that measure- 
ments of the cosmic microwave background radiation can 
be used to derive an upper limit on the gravitational- 
wave spectral energy density parameter, na mely , ficw < 
1.4 X IQ-io for ujo = 1-4 x lO'^^ rad/s (Hf. Some 
gravitational-wave spectra shown in this paper do not 
satisfy this bound. Taking into account that the inflaton 



mass 7710 determines the overall vertical displacement of 
the gravitational-wave spectrum, one just needs to con- 
sider lower values of in order to make these spectra 
compatible with the above mentioned upper limit. 



IV. CONCLUSIONS 

In this work we have investigated the generation 
of gravitational waves within loop quantum cosmology 
models. For such models, the evolution of the universe 
is naturally divided in two stages, corresponding to the 
semiclassical and the classical regimes. In the former, 
loop quantum effects introduce modifications to the dy- 
namical equations describing the evolution of the early 
universe, while in the latter the evolution proceeds ac- 
cording to the usual general relativity equations. The 
transition between the two regimes takes place at very 
early times. 

For the semiclassical regime we have assumed that the 
corrections to the dynamical equations are of the inverse- 
volume type, leaving for a future investigation the holon- 
omy corrections. Inverse-volume corrections involve two 
ambiguity parameters, j and £, which we have assumed 
to be free parameters. 

For the classical regime, we have assumed that the evo- 
lution of the universe proceeds according to the usual 
standard inflationary model, i.e., a inflationary epoch (of 
the chaotic type) is followed by reheating and then, suc- 
cessively, by radiation-dominated, matter-dominated and 
dark energy-dominated periods of evolution. 

We have also assumed that, initially, the inflaton field 
(/) is located near the origin of the chaotic-type poten- 
tial, taking values much smaller than the Planck mass. 
This choice of initial conditions is admissible, since loop 
quantum corrections to the dynamical equations guar- 
antee that the scalar field (f) is pushed up the potential, 
reaching maximum values of the order of the Planck mass 
at the beginning of the inflationary period (see Fig. [1} , 
which is required to obtain enough inflation. 

The assumption that ^ mp and ^ mp, to- 
gether with the requirement that initially the scalar field 
satisfies marginally the uncertainty principle, imposes a 
lower bound on the value of the ambiguity parameter j, 
namely, j > 100. Another constraint on the values of j 
and £ comes from the requirement that the scale factor 
grows at least 60 e-folds during the inflationary period. 
Our numerical calculations show that this condition is 
satisfied just for the values of the ambiguity parameters 
j and £ corresponding to the shaded region of Fig. 

Loop quantum effects introduce modifications not only 
to the dynamical equations describing the evolution of 
the universe, but also to the equation of tensor modes. 
However, in this paper, we have taken into account only 
the modifications to the dynamical equations, thus sim- 
plifying considerably the calculations without loosing the 
main features of the gravitational- wave energy spectrum. 

To calculate the full gravitational-wave energy spec- 
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trum we have used the method of continuous BogoHubov 
coefficients. Our analysis shows that, for certain condi- 
tions, loop quantum effects leave a clear signature on the 
spectrum, namely, an over-production of low-frequency 
gravitational waves. This signature is present on the 
gravitational-wave spectrum only if the growth of the 
scale factor during the inflationary period docs not ex- 
ceed significantly the minimum growth required in stan- 
dard inflationary cosmology, namely, 60 e-folds. If the 
scale factor grows much more than this value, gravita- 
tional waves generated prior to the inflationary period 
have, today, a wavelength much bigger than the Hub- 
ble horizon, leaving no imprint on the gravitational- wave 
spectrum. On the other hand, if the growth of the scale 
factor during the inflationary period is smaller than 60 
e-folds, then the imprint of loop quantum cosmology is 
clearly seen on the gravitational-wave spectrum, at fre- 
quencies which increase with decreasing number of e-folds 
of expansion during inflation. For example, in the case 
I = 0.6 and j = 100, for which about 28 e-folds of expan- 
sion are obtained during inflation, loop quantum effects 
leave their signature at LISA frequency band. Taken 
into account the above comments, we conclude that the 
values of the ambiguity parameters j and i for which a 
signature of loop quantum cosmology shows up on the 
gravitational-wave spectrum are those corresponding to 



a narrow band around the thick line of Fig. [2l i.e., these 
values for which the scale factor grows about 60 e-folds 
during the inflationary period. 

Our results shows that, contrary to what is usu- 
ally assumed, inflation does not necessarily erase all 
the information on physical features present in the pre- 
inflationary era. Indeed, as we have shown, within loop 
quantum cosmology physical processes taking place in 
the very early universe, prior to the inflationary period, 
leave their imprint on the spectrum of the gravitational 
waves at very low frequencies, corresponding to wave- 
lengths of the order of the Hubble distance. Despite the 
fact that the gravitational-wave spectral energy density 
parameter ^2gw for such frequencies may be quite high, a 
direct detection is not possible. Nevertheless, these grav- 
itational waves may have left their imprint on the cosmic 
microwave background radiation and on large scale struc- 
tures, in which case they will allow us to test present-day 
theories about the quantum origin of the universe. 
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